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The charge and spin response of a spin-polarized electron gas is investigated including terms beyond 
the random phase approximation. We evaluate the charge response, the longitudinal and transverse 
spin response, and the mixed spin-charge response self-consistently in terms of the susceptibility 
functions of a non-interacting system. Exchange-correlation effects between electrons of spin a 
and cr are included following Kukkonen and Overhauser, by using spin-polarization dependent 
generalized Hubbard local field factors Ga-^ and G^^. The general condition for charge-density and 
spin-density- wave excitations of the system is discussed. 
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Response functions relate the induced charge and spin 
densities to the strength of an external disturbance and 
play an important role in the understanding of many- 
body systems. The spin-polarized electron gas (SPEG) 
is an n-electron system with electrons of spin a and 
electrons of spin a embedded in a uniform positive charge 
background. (The volume of the system is taken to be 
unity in this work.) Previous investigations of the re- 
sponse of the spin polarized electron system were limited 
in scope. Some focussed on the paramagnetic responsetl 
or on the charge and spin-density fluctuations of a fer- 
romagnetic electron gas within the Hartrce-Fock(HF) 
approximation^. Others used the random phase approx- 
imation (RPA)0, local spin-density functional theoryo, 
or were limited to the infinitesimally polarized electron 
liquidu. However, the role of correlations beyond the 
RPA in the charge-spin response has never been exam- 
ined explicitly for the case of arbitrary spin polarization 
< Id < 1, where C = {n^ - n^)/ina- + n^). 

The purpose of this paper is to present a treatment 
of charge and spin response in a unified way. The self- 
consistent effective potential experienced by an electron 
of spin s is expressed in terms of the charge density fluc- 
tuation Sn and the spin density fluctuation Sfn. The 
exchange-correlation interactions between electrons of 
the same spin (ss) or of opposite spins (ss) are included 
by employing Hubbard-type spin-dependent local-field 
factors Gs^ and Gg^ . The self-consistent linear response 
method of Kukkonen and OverhauserS is extended to a 
SPEG by generalizing the local field factor. The charge 
and spin response to an arbitrary electric and magnetic 
disturbance is derived and compared with the existing 
theories. 

We consider an electron gas in the presence of a uni- 
form positive charge background. The imbalance in the 
populations of up and down spins forming a system of 
SPEG is caused by an effective dc magnetic field whose 
origin needs not be specified in detail. Any degree of spin 
polarization C can be obtained by adjusting the value of 
the effective magnetic field B. We assume the SPEG is 



disturbed by an infinitesimal external electric potential 
^^ext^^ and magnetic field bo{r,t). In response to these 
external electric and magnetic disturbances, charge and 
spin fluctuations are set up in the system, and the Hamil- 
tonian for an electron with spin s can be approximated 
as 



H = Hn 



Hi' 



(1) 



where Hq is the Hamiltonian of a single quasiparticle of 
the SPEG in the absence of the external disturbance. 
Hi'^ is the spin-dependent self-consistent effective per- 
turbation. The eigenstates and eigenvalues of Hq are 
given by |fc, (J > and ecr(fc). In this work we assume the 
spin-splitting is much greater than the Landau level split- 
ting and ignore any degree of orbital quantization. Since 
the most general disturbance can be decomposed into its 
Fourier components, we choose the disturbances Uo"^^*, 
and Hi'^ to vary as ~ ^i^t-iq-r ^ The self-consistent mag- 
netic disturbance b is the sum of &o and 47rm, where m 
is the induced magnetization. The Fourier component of 
the most general Hi'^{f,t) can be written asQ 



las-b + vo"'-'* +v{q)[5n (1 



G, 



(2) 



In Eq. (g), for the sake of brevity, the q and uj depen- 
dence of the local flelds, fluctuations, and disturbances 
has not been displayed. The parameter 70 is given by 
70 = 53*/^ -B with g* , hb, s, and v{q) being the effective 
g-factor, the Bohr magneton, Pauli spin operator, and 
the Fourier transform of the bare Coulombic potential, 
respectively. Equation (^) is the generalization of the ef- 
fective interaction Hamiltonian of the SPEG in the pres- 
ence of inflnitesimal magnetic and electric disturbances. 
The local fields Gs^ and Gs~ are responsible for charge- 
d spin-induced correlation effects on an electron of spin 



Gs — Gc 



±G, 



G.± = Gs 



±Gs 



(3) 
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where Gss^'^ {Gss^'^) and Gss'^ [Gss'^) account for the 
paraUel-spin exchange-correlation and the antiparallel- 
spin correlation effects in linear response theory. 

The charge- and spin-density fluctuations Sn{q,Lu), 
S'mi{q,u}) {i = z,+, and — ) are given, in terms of the 
self-consistent effective perturbation w), using the 

equation-of-motion of the density matrixB. 

By taking the matrix element of the effective pertur- 
bation, Eq. (H), with respect to eigenstates s >, then 
Fourier transforming the resulting expressions and com- 
bining the results with the deflnitions of charge- and 
spin-density fluctuations, we obtain the coupled equa- 
tions for < si\Hi''(q,uj)\s2 > in terms of the external 
charge and spin disturbances. We then solve the coupled 
equations for < si|iJi'*|s2 > and substitute the corre- 
sponding matrix elements back into the expressions of 
the fluctuations. The charge-density fluctuation lon- 
gitudinal spin-density fluctuation Sniz, and transverse 
spin-density fluctuations Sm+ and Srri- can then be ex- 
pressed in terms of a susceptibility matrix x as 



—eSn 



boz 



and 



X-""" / Uo" 



(4) 



(5) 



X-'"™(g,^) = — — -f 

1 + ^(vGg- 



{vGg - 47r7o) Use 



(11) 



In Eqs. (|) to jnl), Ilsisa" is the Lindhard-type electric 
(si = S2) or spin (si 7^ S2) susceptibility. 

n.,.,°(g,c.)=V nois..Ak + ql)-Me.Ak)) ^^^^ 
Esjk + q) - es2ik) - huj + 17] 



where no(£s(fc)) denotes the equilibrium distribution of 
quasiparticles having spin s. The D is given by 

D{q,u;) = i[l - 2vU,,°{l - - Gg,^)] 

[1 + 2vIlgg°{Gsa''' - G^sl - 87r7on^5°] 

[1 + 2wn,,°(C"= - Gsa') - 87r7on,.0]. (13) 

The various terms containing factors proportional to tt^q 
in the expressions of susceptibilities have their origin in 
the use of the self-consistent magnetic disturbance. If we 
neglect the induced magnetization m, those terms dis- 
appear from the expressions for various susceptibilities. 
The spin-polarization dependent Fermi wave number of 
the majority(minority) electrons with spin is given 

by 



Here 60 — box ± iboy, and t/Jo'^^* denotes the external 
electric potential corresponding to the external distur- 
bance ?;o'^'^*(= —eifo''^*)- The various susceptibilities are 
written as 



2t;n,,on^/(G^- + G,- 



X^^{q,u;) = -^[U,,'-U,,^ 
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(6) 



+2t;n,,"n^^"(G^--G,-)], (7) 



2,;n<,<,"n^/(G^+-G,+)], (8) 



2t;n„,"n^5"(2-G^+-G,+)], (9) 



l-Ki(z;G,- -4^70) n^ff"' 



(10) 



and 



Af„(s> = fcFo(l±C) 



1/3 



(14) 



where hp^^ is the Fermi wave number for the unpolarized 
case. The expression of ^ss^ is obtained by replacing 
the quantities of spin indices u in the expression of ^^J^ 
by that of a. The ^as^ and ^sJ^ appearing in Eqs. ( p^ ) 
and (^ij) are the susceptibility functions of the spin-flip 
processes. 

In the absence of the perturbation, the charge den- 
sities associated with the majority- and minority-spin 
electrons are spatially uniform but unequal. Hence they 
have only a non- vanishing g = Fourier component. 
However, in the SPEG a spatially varying electric or 
magnetic disturbance with flnite wave number q induces 
electron density fluctuation of each spin, bua and (Jn^, 
and hence a finite spin density fluctuation 8fh. The x*^^, 
and x±""" are the ordinary charge-, longitudinal 
spin-, and transverse spin-susceptibilities, respectively. 
The off-diagonal mixed susceptibilities x'^™(x™'^) corre- 
spond to the charge density response to a magnetic dis- 
turbance (the longitudinal spin density response to an 
electric disturbance). The susceptibilities given by Eqs. 
(^ to (11) reduce to appropriate forms for the unpolar- 
ized or infinitesimally polarized limits. If we set all the 
local fields Gss' = 0, we obtain the RPA susceptibilities 
of the spin-polarized system. When C = 0, Hcrcr° — ^ss^ 
and Ga^ ~ Gs^ . In this case, y^^ and Xlh"™' respec- 
tively, reduce to the well known expressionsQ, 
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l-t;no°(g,^)(l-G, 



(15) 



and 



" ' l + no°(g,c.)(z;a--4^7o)' ^ ^ 

where Uo° = ^aa" + ^ss° = 2U^^°. Tie mixed re- 
sponses vanish in the unpolarized systenm. But, they 
become finite and equal to eafih other in the case of in- 
finitesimally polarized systemQ, in which H^a'^ ^ ^ss^ 
but Gcr^ = Ga^. However, since it is not necessary 
that GcTCT*^ — Gga'^ in the system with finite spin polar- 
ization, we conjecture from Eqs. (Q) and (H) that, for 
the most general case, the charge-spin cross susceptibil- 
ities x*^™ a-nd x™*^ could be different. The inequality of 
the cross-correlation local fields Go-a-^, Gaa'^ is expected 
from the fact that the density n\j'^ (r) , associated with the 
exchange-correlation hole around a given electron with 
spin cr located at the origin, and with n^^^ (r), that around 
an electron with spin ct, are given in terms of the corre- 
sponding pair correlation functions ^o-o-' {'>') by 



(r) 



(17) 



nf\r) = n^[l - g^^(r)] + n,[l - g^^ir)]. (18) 

Equations (0) and (|l|) show that the ni'*' (r) ^ nf^ (r) 
in the SPEC Within the HF or RPA-type approxima- 
tions. Gas'' = = Gga", heucc, wfiJjiave that x*"" = X™*" 
even in the spin polarised systemHa'EI. In the local spin- 
density approximationQ, Gunnarsson and Lundqvist ob- 
served that x*^™ = X""^ by keeping only the diago- 
nal elements of the matrix C_ in their expressions for 
the offpdiagonal charge-spin susceptibilities (Eq. (69) 
in Ref.cl). The matrix elements Cs.s' in Ref.cl are di- 
rectly related with the local fields Gs{s)^ in the present 
work. In the same context, Eqs. (10) and (11) suggest 
that, in general, the transverse spin response functions 
X-i-'"™(9,w) and X-'""'(9ji£) could be different in the 
SPEC Equations (||) to ( [III ) can be considered as def- 
initions of the wave number- and frequency-dependent 
local fields Go-^ and Gg^^, in terms of the corresponding 
fluctuations, in the SPEC Within the HF approxima- 
tion the local fields satisfy the relation G^'^ = Ga~ and, 
hence, the mixed charge-spin response functions become 

equal(x*'" = x""")- 

Because the divergences in these response functions 
give the collective modes in the system, various suscepti- 
bility functions obtained here can be used to investigate 
the collective modes such as charge-density and spin- 
density wave excitations in the SPEC The coupling of 
charge-density waves and spin-density waves is expected 
in the SPEG, and the conditions for the spin-flip trans- 
verse modes are written, from Eqs. (10) and ([TTI), by 
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[vGa - 47r7o) H 



0, 



(19) 



l + ^{vGs- -At:^o) Hff^° = 0. 



(20) 



For example, the divergence of leads us, in a long 

wavelength limit, to a mode 

nC 

hu = 2joB + ^(47170 - vGg-) + a(C) <z^ (21) 

where the coefficient a depends, in general, on the degree 
of spin polarization of the system and reduces to a(l) = 
^(^^ bm^okp ) ^'^^ complete spin polarization 

C = 1. The second term of Eq. (|2^) disappears as the spin 
polarization of the system vanishes. On the other hand, 
the general expression for the dispersion relation of the 
coupled longitudinal modes is given by the zeros of the 
D{q,u;) defined by Eq. 

= [1 - t;H,,"(l ~ 2C^^)][1 - vUgg'^il - 2Ggsn] 



- 47r7o(no 



H, 



87r7own,,°Hss"(2 - G,+ - Gg- 



- i;2n,,"Hffs"(l - 2G.s")(l - ^Gsa"). 
In RPA, the above expression reduces to 

i_(i, + 4^^o)(n,,Vngg°) 

+ 16^i;7oH^,°n5g° = 0, 



(22) 



(23) 



where Hcrcr° 7^ ^ss'^ ■ One can expect, from Eq. ( [2^ ) 
coupled modes of charge-density and spin-density wave 
excitations, with a long wavelength limit dispersion rela- 
tion given by 

LJ+{q) =UpL + —- [ (1 - — ^)-2 



rig 

qTFi, 



uj^ (q) = CI 



PL\ 



An„ng 70 



9 , Tl(j Tl(j 



5n^ qTpJ qTFs 



q. 

(25) 



Here flpL is the plasma frequency corresponding to the 
total electron density n = + ng and qTF^f^^j is the 
Thomas-Fermi wavenumbert2l of the majority(minority) 
electrons of spin cr(a). The terms involving 70 have their 
origin in the self-consistent magnetic responses and, es- 
pecially, those terms containing Anang result from the 
coupling of the electric and magnetic responses in Eq. 
(p3|). The terms of rio-^ and Ug"^ are due to the contri- 
bution to the noninteracting susceptibility of the extra 
kinetic energy in the SPEG Fermi sea. The minus sign 
in Eq. ( p5| ) indicates that, in the long wavelength region, 
the slope of the dispersion of the mode omega- is sup- 
pressed in the presence of the coupling between the oscil- 
lations of charge density and spin density of the SPEG. 
When C = 0, Eq. (22) reduces as 
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[1 - vIlo°{l - G,+)] [1 + n„°{vG,- - 4^70)] = 0, (26) 



which is the product of the conditions of self-sustaining 
oscillations of charge and spin densities in a spin- 
unpolarized system given by Eqs. ( |l5|) and (^6|). For 
the case = 0, Eq. (||) reduces to [1 - vno"][l - 

47r7ono°] = 0, which is the RPA result of a unpo- 
larized system. The first factor leads us to the well- 
known charge density wave excitation due to Coulomb 

interactionll3; to — ^pl{^ + nigTT^)- ^"^^ other 
hand, the second factor gives us the spin density wave 
excitation, in response to the self-consistent magnetic 
disturbance, of the spin-unpolarized Fermi sea; to — 

f^PL\/5'(l + §7^977^)^^^ ■ "^^"^ '^^^^ °^ complete 
spin-polarization C = 1; the frequency and wave number 
dependence of various longitudinal susceptibilities given 
by Eqs. to (||) becomes the same and the condition 
for the longitudinal collective modes is given by 

l-n,/[i;(l-2G,-^)+47r7o] -0. (27) 

If we set Ilaa" + n^ffO but G^^ = Gs^ , Eq. (|2|) be- 
comes 

= 1- i;(n,," + ns5°)(i - G,+) + [n,,° + n^^" 

- 4i;n,,''n5s"(l - G,+)](«G,- - 4^70). (28) 

Taking the external magnetic disturbance &o as our ef- 
fective magnetic disturbance in Eq. (||), instead of the 
self-consistent field &, makes the factor 47r7o on the right 
hand side of Eq. ( p8| ) disappears, and the expression re- 
duces to the result of an infinitesimally spin polarized 
systemEl. 

In summary, a unified treatment of the response of 
the spin polarized electron gas is presented in this paper 
and general expressions for various susceptibility func- 
tions are derived. The present results reproduce exactly 
the known results for several simple situations. Spin- 
polarization dependence of the HE local fields is displayed 
explicitly. We believe that our results could be useful in 
understanding electric, magnetic, and optical properties 
of a number of spin polarized systems like ferromagnetics 
and diluted magnetic semiconductors. 
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